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ABSTRACT 



The main goal of this work is to calculate the contributions to the cosmological recombination spectrum due to bound-bound transi- 
tions of helium. We show that due to the presence of helium in the early Universe unique features appear in the total cosmological 
recombination spectrum. These may provide a unique observational possibility to determine the relative abundance of primordial 
helium, well before the formation of first stars. We include the effect of the tiny fraction of neutral hydrogen atoms on the dynamics of 
He II — > He I recombination at redshifts z ~ 2500. As discussed recently, this process significantly accelerates He ll ^ He I recombi- 
nation, resulting in rather narrow and distinct features in the associated recombination spectrum. In addition this process induces some 
emission within the hydrogen Lyman-cc line, before the actual epoch of hydrogen recombination round z ~ 1 100- 1500. We also show 
that some of the fine structure transitions of neutral helium appear in absorption, again leaving unique traces in the Cosmic Microwave 
Background blackbody spectrum, which may allow to confirm our understanding of the early Universe and detailed atomic physics. 

Key words, atomic processes - cosmic microwave background - cosmology: theory ~ early Universe 



1. Introduction 

The recombination of helium practically does not influence 
the Cosmic Microwave Background (CMB) an gular fluctua- 
tions, as measured with great success by Wmap dBennett et alj 
l2003h . since it occurred well before the Thomson visibility func- 
tion d efined by hydrogen recombination (Sunyaev & Zeldovich 
Il970l) reaches its maximum. However, similar to the release 



of photons during the ep och of cosmo l ogical hydrogen re- 
combination ( Rubino-Mart in et all 120061: IChluba et alT 120071: 
IChluba & Sunvaevil200 6a). one does expect some emission of 
photons by helium, and the main goal of this paper is to calcu- 
late the contributions to the cosmological recombination spec- 
trum due to bound-bound transitions of helium. 

In our recent papers we computed the detailed cosmo- 
logical recombination spectrum of hydrogen resulting from 
bound-bound (Rubino-Martin et al. 2006; Chluba et al.' | 2007| ) 
and bound-free (Chluba & Sunvaev 2006a) transitions between 
atomic levels, including up to 100 shells, also taking into account 
the evolution of individual energetically degenerate angular 
momentum sub-states. We follo wed the ideas a nd suggestions 
of earlier i nvestigations (Z eldovich et al.l Il968t [Peebles 1968; 
DubrovichI 1975; Bernshte in et alj 1977 : Beigman & Sunvaev 



19781; iRvbicki&deir Anto nio" 1993^, 'Dubrovich & Stolyarov 



19951; iBurginI 120031; iDubrovich & Shakhvorostova 2004 



Kholupenko et alj 120051 IWong et all 120061) . Observations of 



these recombinational fines migfit provide an additional un- 
biased way to dire ctly determi n e tfie baryo n density of tfie 
Universe (e.g. see iDubro vicn ('1975') and 'Bernshte in et alj 
(Il977h . or more recently Kholupenko et al. (2005)) and 
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IChluba & Sunyaev! (l2007bl) ) and to obtain some additional 
information about the other key cosmological parameters, 
facing different degeneracies and observational challenges. 

Obviously, direct evidence for the emi ssion of extra ~ 5 
photons per recombining hydrogen atom dChluba & SunvaevI 
2006a) will be an unique proof for the completeness of our un- 
derstanding of the processes occurring at redshifts z ~ 1400, 
i.e. before the CMB angular fluctuations were actually formed. 
From this point of view an observation of lines emitted during 
He III -» He ii close to z ~ 6000, and He ii -» He i around 
z ~ 2500 will be an even more impressive confirmation of 
the predictions within the standard hot big bang model of the 
Universe, realising that nowadays exact computations using the 
full strength of atomic physics, kinetics and radiative transfer in 
principle should allow a prediction of the cosmological recom- 
bination spectrum from both epochs with very high precision. 

The first attempt to estimate the emission arisin g from 
helium recombination was made by Dubrovich & Sto lvarovl 
(1 19971) . However, only now detailed numerical computations 
are becoming feasible, also due to the fact that atomic physi- 
cists be gan to publish accurate and user-friendly transition 
rates (IDrake & MortonI 120071: iBeigman & Vainshtein ' 2007) for 
neutral helium, including singlet-triplet transitions, which very 
strongly influence the recombination of helium. 

According to the computations of nuclear reactions in the 
early Universe (Olive & Steigman 1995; Cvburt 2004), the abun- 
dance of helium is close to 8% percent of the number of hydro- 
gen atoms, so naively only small additional distortions of the 
CMB blackbody spectrum due to helium recombination are ex- 
pected. However, for helium there are two epochs of recombina- 
tion, a fact that at least doubles the possible amount of additional 
photons. Furthermore, He ill — > He ii recombination is very fast, 
in particular because there is a large quasi-constant amount of 
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free e lectrons belonging to hydrogen dPubrovich & Stolvarovl 
Il997l) . This implies that photons are emitted in a much shorter 
period, so that more narrow features are produceqj. It is also very 
impressive that the He in —> He ii recombination lines practi- 
cally coincide and therefore amplify the corresponding hydrogen 
line (see Fig.[T]i- This is because the difference in the redshifts of 
the two recombinations is close ~ 4.3, when on the other hand 
the energy of similar transitions scales as Z^ = 4 for He ii, such 
that the two effects practically compensate eachother. 

The spectral distortion due to He ii ^ He i recombination 
should have a completely different character First, for small n 
neutral helium has a much more complicated spectrum than hy- 
drogenic atoms (e.g. highly probable fine-structure transitions). 
In addition, the ratio of the energies for the second and first shell 
is ~ 2.1 times higher than for hydrogenic atoms, while the en- 
ergies of the highly excited levels are very close to hydrogenic. 
Since the transitions from the second to the first shell are control- 
ling helium recombination, this leads to the situation that even 
for transitions among highly excited levels the corresponding 
An = 1 -lines do not coincide with those emitted during hydrogen 
or He III -^ He ii recombination. 

Also it will be shown below (Sect. |5]l, that in the recombi- 
national spectrum some fine-structure lines become very bright 
and that two of them are actually appearing in absorption. These 
features lead to additional non-uniformities in the spectral vari- 
ability structure of the total CMB spectral distortion from re- 
combination, where some of the maxima are amplified and oth- 
ers are diminished. This may open an unique possibility to sepa- 
rate the contributions of helium and hydrogen, thereby allowing 
to measure the pre-stellar abundance of helium in the Universe. 
Until now not even one direct method for such a measurement is 
known. 

For the computations of the recombinational helium spec- 
trum we are crucially dependent in the recombination history 
of helium and additional processes that affect the standard pic- 
ture strongly. In this context, probably the most important phys- 
ical mechanism is connected with the continuum absorption 
of the permitted 584 A and intercombinational 591 A line by 
a very small amount of neutral hydrogen present in ioniza- 
tiona l equilibrium during the time of He il -^ He I recombina- 



tionjHuetal.11199 5': 'Switzer & Hirata"2007a'; 'Kholupenko et al.' 
120071) . ISwitzer & Hirata (2007a) and Kholupenko et al. (2007) 
recently made the first detailed analysis of this problem, and in- 
cluded it for the computations of the He ii — » He i recombination 
history, showing that the recombination of neutral helium is sig- 
nificantly faster. Here we reanalyse this process, and discuss in 
detail some physical aspects of the escape problem in the afore- 
mentioned lines. 

We first consider two "extreme" cases for the escape prob- 
lem, which can be treated analytically: (i) where line scattering 
leads to complete redistribution of photons over the line profile, 
and (ii) where there is no redistribution Moreover, we develop 
an useful ID integral approximation for the escape probability 
which permit us to treat any of these two cases without increas- 
ing the computation time significantly. Our final results for the 
escape probability in the more realistic case of partial redistribu- 
tion (or equivalently for coherent scattering in the rest frame of 
the atom) are based on detailed numerical computations which 



' As we will demonstrate here, even the scattering of photons by free 
electrons cannot change this conclusion (see Sect. 15. 5t . 

^ In this case line scattering is totally coherent in the lab frame. 
During recombination this is a very good approximation in the very 
distant wings of the line. 



will be presented in a separate paper (Chluba et al. 2007, in 
preparation). One can then obtain a sufficiently accurate descrip- 
tion of the real dynamics hy fudging the escape probability using 
the "no redistribution" case mentioned above with a certain func- 
tion that can be obtained by comparison with the full numerical 
computations. Our final results for the escape probability are in 
very good agreement with those obtained by ISwitzer & Hiratal 
(2007a) (see discussion in Sect.[3]l. 

Interestingly, the hydrogen continuum process leaves addi- 
tional distinct trace in the cosmological recombination spec- 
trum, because the continuum absorption of the He: photons 
leads to significant early emission in the H i Ly-a transition at 
V ~ 1300 GHz, well separated from the Ly-c line originating 
during hydrogen recombination at v > 1500 GHz, and contain- 
ing about 7% of all photons that were released in the hydrogen 
2p-ls transition. The amplitude and width of this feature is com- 
pletely determined by the conditions under which the above pro- 
cess occurs. 

The main result of this paper are the bound-bound spectra 
of Hen and Hei from the epoch of cosmological recombina- 
tion (see Fig. [TJ. The strongest additions to the cosmological 
hydrogen recombination spectrum due to the presence of helium 
lines reach values up to 30-40% in several frequency bands. This 
strongly exceeds (roughly by a factor of four) the relative abun- 
dance ratio of helium to hydrogen, raising hopes that these dis- 
tortions will be found once the recombinational lines will be- 
come observable. 

It is important to note that for the computations in this paper 
(see Fig.[T]i, we do not include the impact of feedback processes 
on the computed recombinational lines. Among all the possible 
feedback mechanisms, the most relevant for the recombinational 
spectrum is the pure continuum absorption (far away from the 
resonances) of the remaining He 1 2^^Pi - 1 ' So intercombination- 
line and He 1 2'Pi - 1 'So photons. Due to this process, these pho- 
tons will be finally absorbed, and the corresponding features on 
the final spectrum will disappear, producing additional photons 
that will emerge mainly through the Lya line in the hydrogen 
spectrum. 



2. Basic equations. The Helium atom 

A description of the basic formalism and equations to calcu- 
late the time-evolution of the populations for different atomic 
species (hydrogen or helium) within a multi-level code during 
the epoch of co smological recombin ation (800 < z < 7000) 
can be found in ISeager et al.l (l2000l) . In this paper, we follow 
the same approach and notation that was used in our previous 
work s for the computation of the hydrogen recombin ation spec- 
trum (iRubino-Martin et al.ll2006tlChluba et cd]|2007h . The codes 
used for the computations presented here were obtained as an 
extension of the existing ones, by including the equations for the 
population of the He i levels. As in our previous works, we devel- 
oped two independent implementations in order to double-check 
all our results. 

For all the results presented in this paper we use the same 
values of the cosmological parameters which w ere adopted in 
our previous works, namely dBennett et al.ll2003b : Ob = 0.0444, 
Qt^, = 1, Q™ = 0.2678, Qa = 0.7322, 7p = 0.24 and /i = 0.71. 

2.1. He I model atom 

In our computations we follow in detail the evolution of the level 
populations within neutral helium, including up to Wmax - 30 
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Fig. 1. Full helium and hydrogen (bound-bound) recombination spectra. The following cases are shown: (a) the He ii —> He i 
recombination spectrum (black solid line), which has been obtained including up to Wmax - 30 shells, and considering all the J- 
resolved transitions up to n = 10. In this case, there are two negative features, which are shown (in absolute value) as dotted lines; 
(b) the Hem —^ Heii recombination spectrum (red solid line), where we include «niax = 100 shells, resolving all the angular 
momentum sub-levels and including the effect of Doppler bro adening due to scattering off free electrons; (c) the H i recombination 
spectrum, where we plot the result from' Chluba et al.f ( l2007h up to n„ax = 100. The Hi Lyman-a line arising in the epoch of Hei 
recombination is also added to the hydrogen spectrum (see the feature around v = 1300 GHz). In all three cases, the two-photon 
decay continuum of the n - 2 shell was also incorporated. Feedback processes for the He ii and He i recombinations are not taken 
into account. Blue line shows the total recombination spectrum. 



shells. For all levels, we distinguish between "singlet" (S - 0) 
and "triplet" (S = 1) states. Up to n = 10, we follow sepa- 
rately all levels with different total angular momentum J. This 
permits us to investigate in detail the fine structure lines ap- 
pearing from cosmological recombination. Above n - 10, we 
do not resolve in J quantum number, and only LS-coupling is 
considered. Each individual level is quoted using the standard 
"term symbols" as n^^^^Lj, and the spectroscopic notation is 
used. When considering J-resolved levels, the degeneracy fac- 
tor is given by gi = (27 -i- 1), while in opposite case we would 
havegi^(2S + 1)(2L+ 1). 

To completely define our model atom, we need to specify 
for each level / = {n,L,S,J], the energy, £,, and photoioniza- 
tion cross-section, cr,c(v), as a function of frequency, which is 
important in order to take into account the effect of stimulated 
recombination to high levels. Finally, a table with the Einstein 
coefficients, Aj^j, and the corresponding wavelengths for all the 
allowed transitions has to be given. 



2.1.1. Energies 

The energies for t he diff erent levels up to « = 10 are taken from 
iDrake & MortonI (|2007[) . However, this table is not absolutely 
complete, and some of the high L sub-states for outer shells are 



missing. In order to fill this table up to n^ax - 30, we proceed 
as follows. We use the formulae for quantum defects (see lDrakg 
1996, Chap. 1 1) to compute the energies of all terms with L < 6 
and n > 10. For all other levels, we adopt hydrogenic values 
for the energies, using -Rii/rP', with Ru «; 13.6 eV. Note that 
in this last case, the energy levels will be degenerate in L and 
J. However, this approximation is known to produce very good 
results (Beigman & Vainshtein 2007). A summary of the final 
energies adopted for each particular level in our He i model atom 
is shown in Figure |2l both for the singlet and triplet states. 

2.1.2. Photoionization cross-sections 

For n < 10, the photoioniz ation cross-section s, cridv), are 
taken from TOPbas^l database dCunto et al.lll99 3'). We note that 
this database does not contain 7-resolved information and only 
cross-sections for L < 3 can be found. To obtain /-resolved 
cross-section we assume that the cross-section for each sub-level 
is identical, i.e. 

There are two important issues that we would like to stress. 
First of all, there are large gaps in the tables from TOPbase. 



^ At |http://vizier.u-strasbg.fr/topbase/topbase.html[ 
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Up to n = 5, we co mpute d the missing cross-section s usin: 
the expressions from Smits (1996) and Benjamin et alj (1199' 
for the spontaneous photorecombination rates to infer the pho- 
toionization rates Rjc- With this procedure it is not possible to 
include the effect of stimulated recombination self-consistently. 
We estimated the errors due to this effect using the replacement 
Rjc — > Ric X [1 + rtbhivc)], where nbb(i'c) is the blackbody pho- 
ton occupation number at the ionization threshold of the level, 
and found chang es of the order of 10 % - 20%. For all other 
levels we follow iBauman et alJ (l2005h and adopt re-scaled hy- 
drogenic cross-section. This is done usi ng our computations of 
Cic(y) for the hydrogen atom (based on Karzas & Latter! 1 19611; 
I S tore V & Hummer 1991) , and shifting the threshold frequency 
accordingly. 

Secondly, we would like to point out that the cross-sections 
provided in TOPbase are sparsely sampled. For example the 
power-law behaviour up to twice the threshold frequency is usu- 
ally given by ~ 10 points. Furthermore, due to auto-ionization 
several resonances exist at large distances above the ionization 
threshold, and many of these extremely narrow features are rep- 
resented by 1 point. Fortunately, because of the exponential cut- 
off from the blackbody spectrum these resonances do not affect 
Ric significantly. Still we estimate the error budget using these 
cross-section to be ~ 10%. 

A summary of the final values adopted for each particular 
level in our He i model atom is also shown in Figure|2] We finally 
note that, in order to speed up our computations, we tabulate the 
photoionization rate during the initialization state of our codes, 
which involves one-dimensional integral over the (blackbody) 
ambient photon field, and we interpolate over this function when 
needed. At every particular redshift, the corresponding photore- 
combination rate is computed using the detailed balance relation, 
which is satisfied with high precision due to the fact that at the 
redshifts of interest, the electron temperature and the radiation 
temperature practically do no differ 

2.1.3. Transition probabilities 

Our basic databa s e for t he transition probabilities is taken from 
iDrake & MortonI (|2007|) . which is practically complete for the 
first 10 shells, and includes 937 transitions between J-resolved 
states. This database does also contains spin-forbidden transi- 
tions (i.e. singlet-triplet and triplet-singlet), which take into ac- 
count the mixing of singlet and triplet wave functions. 

However, there are some transitions missing in these tables, 
which involve lower levels with n > 8 and L > 1 for the sin- 
glet, and n > 7 and L > 6 for the triplet states. These gaps are 
filled using re-scaled hydrogenic values as follows: for a given 
transition {n,L,S,J] —> {n',L',S',J'}, we first obtain the non 
y-resolved transition probability A^l^^i^, scaling by the ratio 
of the transition frequencies to the third power Whenever J- 
resolved information for the level energies is available, we also 
compute the weighted mean transition frequency. However, the 
corresponding corrections are small. To obtain the final estimate 
for the 7-resolved value, we assume that 



(27' + 1) 



nLJ—^n'L'J' 



{lU + \)(2S' + \) 



nL—>n'L' 



(1) 



These expressions are also used to include all transitions involv- 
ing levels with n,n' > 10, and for those between n > 10 and 
n' < 10 states, adopting the corresponding average over J and 
J'. For the case of n^ax = 30, our final model contains 80,297 
bound-bound transitions. 



Finally, we note that once the energies of all levels are ob- 
tained, the wavelengths for all transitions are calculated con- 
sistently using the respective upper (£■„) and lower (£/) energy 
levels as v„/ = (£„ - Ei)/h. This is important in order to guar- 
antee that we recover the correct local thermodynamic equilib- 
rium solution at high reds hifts (see discussion in Sec. 3.2.1 of 
iRubifio-Martin et al.ll2006l) . 



2.1.4. Two-plioton decay and non-dipole transitions 

Apart from the aforementioned transitions, we also include in 
our computations the two photon decays of the 2 'So and 2^Si 
levels. We adopt the values A2iSn-ii So = 51.3s"' andA23s,-iiSo = 
4.09 X lO-"^ S-' (IDrake etal.ll 19691) . 

Our final spectrum also contains the contribution of the 2 'So 
two-photon decay spectrum, which is com puted like for the hy- 
drogen case (see e.g. Eq. 3 in Rubino-Mart m et al.ll2006l) . The fit 
to the profile function (Drake 1986) for this transition is taken 
from Swi tzer & Hirata (2007a). 

We also included some addition al low probability non-dipole 
transitions JLach & Pachuckill200 1'). but in agreement with pre- 
vious studies found them to be negligible. 



2.2. He II model atom 

For singly ionized helium we use hydrogenic formulae (see e.g. 
iRubino-Martin et al...2006) with re- scaled transition frequencies 
(see also ISwitzer & Hiratal l2007al) . The 2s two-photon decay 
profile is modelled using the one for hydrogen, adopting the (re- 
scaled) value of A^J^iJj^ = 526.5 s"'. 



3. Inclusion of the hydrogen continuum opacity 

In order to include the effect of absorption of photons close to 
the optically thick resonant transitions of helium during He ii — > 
He I recombination due to the presence of neutral hydrogen, one 
has to study in detail how the photons escape in the helium 
lines. This problem has been solved by two of us using a dif- 
fusion code, and the results will be presented in a separate paper 
(Chluba & Sunyaev 2008, in preparation). 

For the purposes of this paper, the important conclusion 
is that the results obtained using this diffusi on code are in 
rather good agreement with those presented in Switz er & Hiratal 
(2007a) . showing that for the interaction of photons with the con- 
sidered resonances, the hypothesis of complete redistribution of 
the photons over the Voigt profile, (p{x) (see AppendixlAlfor def- 
initions), is not correct, and may lead to significant differences, 
for the He 1 2'Pi - 1 ' So transition. However, that is not the case 
for the He: 2^Pi - I'So intercombination line, where the real 
dynamics is very close to the full redistribution case. 



3. 1 . The escape probability in the He i lines 

For the computations of the He ii ^ He i spectrum in this paper, 
we make an ansatz about the shape of the escape probability in 
these lines, which is described below. This hypothesis permit us 
to efficiently compute the escape probability in our codes, with- 
out reducing the computational time significantly. This ansatz 
has been tested against the full results (Chluba & Sunyaev 2008, 
in preparation), and is found to produce accurate results for the 
spectrum. To present it, we first describe two particular "limit- 
ing" cases for the escape problem: the complete redistribution 
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Fig. 2. Graphical representation of the values adopted for the energies and photoionization cross-sections of our He i model. For 
each level, two letters are given, which refer to the energy and the photoionization cross sectio n of that leve l, respectively. For the 
energies, the letters refer to: "D": iDrake & MortonI (120071) : "Q": quantum defect expansions Drake (199d); and "H": hydrogenic 
approximation. For the photoionization cross-sections, the letters refer to: "T": TOPbase; "S": iSmits ( 1996) ; Benjamin et al. ( 19991) ; 
and "H": re-scaled hydrogenic values. 



(or incoherent scattering), and the no redistribution case (or co- 
herent scattering in the lab frame). The realistic case will be re- 
ferred as partial redistribution in the line (or coherent scattering 
in the rest frame of the atom). 

For these two particular cases , we foUow the procedure out- 
lined in' Switzer & Hirata (l2007al) . It is based on the assumption 
that within the considered range of frequencies around a given 
resonance a solution of the photon field, including resonant scat- 
tering and hydrogen continuum absorption, can be obtained un- 
der quasi-stationary conditions. Within +1% - 10% of the line 
center (or roughly +10^ - 10'* Doppler width), this approxima- 
tion should be possible!- 

3.1 .1 . Complete redistribution (or incolierent scattering) 



where v(x) = v^'^ +x Avu, N^^ is the number density of hydrogen 
atoms in the Is-state, crfg(v) is the photoionization cross section 
of the hydrogen ground state, H is the Hubble expansion factor, 
and Vq '^ is the transition frequency of the considered helium res- 
onance. The Doppler width, Avp, of the line due to the motion 
of helium atoms is defined in AppendixlAl 

The computational details about the numerical integration of 
Eq. |2] as well as the derivation of a one-dimensional integral 
approximation to the full 2-dimensional integral are discussed in 
Appendix iB] 

3.1 .2. No redistribution (or colierent scattering in tine lab 
frame) 



One finds that in this case, the corresponding correction, Af esc 
to the standard Sobolev escape probability, Ps - [1 - e^'^^]/Ts 
is given by: 



This case corresponds to a situation in which every photon 
coming through the line is emitted again with the same fre- 
quency. During recombination this is a very good approxima- 
tion in the very distant damping wings of the resonance. In prac- 
XM ^oo tise, this case can be treated using the formalism described in 

4>(x)dx I rs0(x')e"^'^*''''*[l -e"^'*'''''']dx'. (2) ISwitzer & Hiratal ([2007a ). For every transition u ^ I. we need 
°° ^ ' to define the following quantity 

Here tl(x, x') = ts [x(x') -xi^)] is the optical depth with respect 
to line scattering off the resonance, where rs is the Sobolev opti- f 

^ (f>(y) dy is the normalized (^(+00) - 1) 
integral over the Voigt-profile. Furthermore we introduced the 
hydrogen continuum optical depth 



R° 



Au^, + R° 



(4) 



Tc(x, X') : 



cN\ 



H Jv(x) y 



where R°'^i is the sum of the rates of all the possible ways of 
leaving the upper level but excluding the considered resonance, 

i.e. ;?r, = R, 



Yjij^iRii^i- In this equation, we have intro- 



(3a) duced the (bound-bound) rates, which are computed as 



cA^f^o-»(v) V 
Hv ^3 



1 



(^r 



(3b) 



R„ 



Ai^uigi/guhbhiyiu), Ei > E„ 



(5) 



'' In addition, the used approximation dyvNy ~ vgdyNy, where vg is 
the transition frequency of the considered resonance and Ny = ly/hv, 
demands that the obtained solution is only considered sufficiently close 
to the line center. For this reason the term connected with emission of 
photon due to the recombination of hydrogen (see definition of Ic in 
ISwitzer & Hiratal jlOOla.) ) should be neglected, since in this process 
practically all photons are emitted very close to the ionization frequency 
v^ «: vo of hydrogen. 



and the photoionization rate, R,i^c. This quantity, fu^i, gives the 
fractional contribution to the overall width of the upper level 
of all possible transitions leaving the upper level except for the 
resonance. In other words, /„^/ represents the branching frac- 
tion for absorption of a line photon to result in incoherent scat- 
tering. During helium recombination /„^; ~ 10"^ for the He: 
2'Pi - I'So transition and close to unity for He: 2^Pi - I'So 
intercombination line. 
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Fig. 3. Contributions to the escape probability of the Hei 
2'Pi - I'So transition (upper panel) and the Hei 2^^Pi - I'So 
intercombination-line (lower panel). Ps denotes the standard 
Sobolev escape probability. The correction to the escape prob- 
ability due to the hydrogen continuum opacity is shown in sev- 
eral approaches: (a) full redistribution case: the full 2D-integral, 
AP2D, as given by Eq. (|2]i, and the ID approximation, APid, ob- 
tained with Eq. ( IB.3l i: (b) no redistribution case: the ID approx- 
imation Pid; and (c) partial redistribution case: the fudged es- 
cape probability based on the ID integration of the no redistribu- 
tion case. For compari son we show the simple analytic approx- 
imations of Kholupenko et al.| ([2007) , and the points extracted 
from Fig. 11 of Switzer & Hirata (2 007al) for the case labelled 
as coherent (which is the equivalent to our partial redistribution 
case) for the upper panel, and the points extracted from Fig. 4 of 
ISwitzer & HirataljioOTbt) . 



Once we have obtained this quantity for the considered tran- 
sition, the corresponding escape probability in the case of fully 
coherent scattering is given by 



fu^lP 



i-(i-/;w)P 



where 



P = ^s(/«^/rs) + APe,sc(/;,^/Ts,Tc) 



(6) 



(7) 



where Ps(fu^iTs) means that the Sobolev escape probabil- 
ity is evaluated at fu^iT^ instead of at ts; and obtaining 
APesc(/«^;Ts,Tc) reduces to the use of equation|2] but evaluat- 
ing it at /a^;Ts instead of at ts, while t^ remains unchanged. 



3.1 .3. Partial redistribution case. Our ansatz 

The detailed treatment of the problem with partial redistribu- 
tion is computationally demanding. Our results are based on a 
diffusion code (Chluba & Sunyaev 2008, in preparation), which 
requires ~ 1 day on a single 3 GHz processor to treat one cos- 
mology. The o ther approach of this pro blem, based on a Monte 
Carlo method dSwitzer & Hiratall2007al) is equally demanding. 

For the computations of this paper, we propose and test an 
ansatz which permit us to compute efficiently the escape proba- 
bility. Our basic assumption is that the ratio of the escape proba- 
bility in the complete problem (partial redistribution case) to the 
escape probability in the problem with no redistribution is a con- 
stant number for a given redshift, or equivalently, it has a very 
small dependence on the recombination history. In that case, we 
can use this function (the ratio of those two cases) to fudge the 
real escape probability in our code in a very fast way. The im- 
portant thing is that we only need to compute a single solution 
of the complete problem in order to tabulate the fudge func- 
tion. Moreover, the reference case (the no redistribution case) 
is fully analytic, and using our ID integral approximation de- 
scribed in Appendix iBl it is obtained very fast. Summarising, 
this scheme permit us to compute the escape probability with 
high accuracy in our codes, without the need of interpolating us- 
ing pre-computed tables. 

In practise, we use the solution for the recombination his- 
tory which was obtained within the no redistribution approxima- 
tion, and we compute the corresponding escape probability for a 
given cosmology. If we make a further iteration, by recomputing 
the new recombination history using the new escape probability, 
we find that the result practically does not change. For the He i 
2^Pi - 1 ' So intercombination-line escape probability is always 
close to the full redistribution case, so for this line we directly 
consider this approximation for the computations. 



3.2. Results for AP^sc 

In this paper, we only consider the corrections to the escape 
probabiUty for the He 1 2'Pi - 1 'So transition and He 1 2^Pi - 1 ' So 
intercombination-line. In principle, all the other n'Pi - 1 'So and 
spin-forbidden transitions are also affected, by the presence of 
neutral hydrogen, but the effect is smaller, and we omit these 
additional corrections for the moment. 

In Fig.|3]we show different contributions to the escape prob- 
ability of the two considered transitions, computed within dif- 
ferent approximations discussed in the last subsection. In our 
computations, for the He 1 2' Pi - I'So transition the effect of hy- 
drogen is starting to become important below z ~ 2400 - 2500, 
whereas in the case of the Hei 2^Pi - I'So intercombination- 
line the escape probability is strongly modified only at z < 
1800 - 1900. One can also clearly see, as illustrated for the 
complete redistribution approach, that in both cases the ID- 
approximation works extremely well at nearly all relevant red- 
shifts. In particular, the differences are small where the devia- 
tions between the inner integrand in Eq. |2] and its analytic ap- 
proximation deduced from Eq. ( |B.3| l are small (see Fig. IB. 11 and 
the discussion in Appendix IbTi. 

For the He 1 2 ' P i - 1 ' So transition, the departure of the escape 
probability with respect to the full redistribution case is very im- 
portant, being at least one order of magnitude different at red- 
shifts below z ~ 2200. Moreover, the full redistribution case be- 
comes important at earlier redshifts (z ~ 2600 - 2700), thus pro- 
ducing a recombination dynamics which would be much closer 
to the Saha solution. In other words, the assumption of complete 



J. A. Rubino-Martin et al.: Cosmological Helium Recombination 



redistxibution significantly overestimates the escape rate of pho- 
tons from the He i 2' Pi - I'So transition, and thus would artifi- 
cially accelerate He ii — > He i recombination. The final (fudged) 
solution is in reality much more close to the "no redistribution" 
case, although the differences with respect to this later case are 
still significant (roughly a factor of 2 at z ~ 1730). Comparing 
our results with other recent computations, we find that our final 
(fudged) solution is very close to the Switzer & Hirata (2007^ 
computation, which was based on a Monte Carlo analysis of the 
escape problem. There are still small differences around red- 
shifts z ~ 2200 - 2400, which could be probably due to the 
fact that we do not include the modified escape for higher lev- 
els. However the formula given in Kholupenko et alj (2007) only 
works at very low redshifts. 

For the He i 2^Pi - 1 'So intercombination-transition, the sit- 
uation is different. The computations, based on the diffusion 
code, show that for this line one can approximate the escape 
probability using complete redistribution at the level of ~ 10%. 
Therefore, for the computations in this paper, we adopt this ap- 
proximation for this particular transition. The lower panel of 
Fig.[3]also shows the comparison between our escape probability 
and thos e obtained in othe r recent publications. The agreement 
with the ISwitzer & Hiratal (12007a) result is again very remark- 
able, except for the redshift region around z ~ 1900. However, 
we have checked that this difference is mainly due to the assump- 
tion of using the full redistribution solution for this line, and that 
this difference implies only small changes in the final recombi- 
nation spectrum. For this transition the actual corrections due to 
electron scattering, which we neglected so far, are larger. Finally 
we note that, although in this case there is an apparent agree- 
ment at low redshifts with the Kholupenko et aL (2(X)7) result, 
their computation corresponds to the quantity Af esc- Thus, when 
adding the contribution of the Sobolev escape, they have a value 
of the probability which exceeds unity. 

3.3. Inclusion into the multi-level code 

In order to account for the effect of the hydrogen continuum 
opacity during He ii ^ He i recombination into our multi-level 
code several changes are necessary. The first and most obvious 
modification is the replacement of the Sobolev escape probabil- 
ity Ps^ Ps + APesc for the He 1 21?, - 1 'So and He 1 2^1 - 1 ' So 
intercombination-transition. Due to the above replacement more 
electrons are reaching the ground state of neutral helium, but no 
additional helium photons are released. Therefore in the compu- 
tation of the helium spectrum the increase in the photons escape 
rate by APgsc should not be included. 

Given the usual net radiative transition rate Ps x ARnis from 
level / to the helium ground state, the increase in the net transi- 
tion rate due to the presence of neutral hydrogen atoms is given 
by AR'*f - APesc X A/?,iis. Since the corresponding photons 
associated with this transition are ionizing hydrogen atoms one 
has to add the rate A7?*f to the electron equation and subtract it 
from the hydrogen Is-equation. Although it is clear that, given 
this small addition of electrons to the continuum, the hydrogen 
ground state population will re-adjust within a very short time, it 
is still possible that the corresponding electrons will reach the 
ground state via various decay channels, including a cascade 
from highly excited levels, which may even end in the 2s level, 
yielding two photons in the two-photon decay transition. Instead 
of assuming that all electrons connected with the increase of the 
net transition rate, A/?*f , are leading t o the emission of a hydro - 
gen Lyman-ff photon only (as done in iKholupenko et alj|2007[) . 
this approach is more consistent. We shall see below that a part 



1.08 



1.06 - 




2000 



2500 



3000 



Fig. 4. The ionization history during the epoch of He ii ^ He i 
recombination for different approaches. 



of the additional electrons indeed take more indirect routes to 
the hydrogen Is-level. 

With these additions to our multi-level code it is possible 
to obtain both the ionization history and the helium and hydro- 
gen recombination spectrum including the effect of the hydrogen 
continuum opacity as outlined in this Section. 

4. The helium recombination history 

The main goal of this paper is to compute the spectral distortions 
resulting from the bound-bound transitions of helium. However, 
since we are discussing several approximations to include the 
hydrogen absorption during the epoch of He ii — > He i recombi- 
nation, we here shortly discuss the corresponding differences in 
the ionization fraction. 

Figure|4]shows our results for the redshift-dependence of the 
free electron fraction x^ - n^jn^ during He ii — > He i recom- 
bination, using the three approximations for the escape prob- 
ability in the Hei lines, as discussed above. Qualitatively, all 
three results (i.e. full redistribution, no redistribution and par- 
tial redistribution) agree with those found in som e earlier studies 
(IKholupenko et al .1120071: 1 Switzer & H irata'200"7a), showing that 
the inclusion of the hydrogen continuum opacity in the compu- 
tation significantly speeds up recombination, making it closer to 
the Saha solution. Our (fudged) solution for the case of partial 
redistrib ution of photons in the r esonance is in good agreement 
with the Switzer & Hiratal (|2007a), except for the small differ- 
ence around z ~ 2200. As pointed out in the last section, these 
are likely due to the fact that we did not include the continuum 
opacity correction for higher transitions. 

It is important to note that the incorrect hypothesis of full 
redistribution in the Hei 2'Pi - I'So -resonance has a strong 
impact on the recombination history. In that case, the effect of 
continuum opacity on the escape probability becomes of impor- 
tance at earlier times, shifting the redshift at which x^ starts to 
depart from the solution without continuum opacity consider- 
ably. In addition, the period during which x^ is very close to 
unity, i.e. just before hydrogen recombination starts, becomes 
considerably longer. Therefore, it is very important for a detailed 
analysis of the recombination history to treat properly the escape 
probability in this line. 

So far we did not consider the effect of fee d back in our 
computations and as shown in ISwitzer & Hiratal (l2007al) one 
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does expect some additional delay of Heii —> Hei recom- 
bination aro und z ~ 24 00. However, looking at Fig. 12 in 
[Switzer & Hiratal (l2007al) . this process is not expected to alter 
the results by more than 10% - 20%. 

Finally, we also mention that for computations of the elec- 
tron fraction during the epoch of He ii — > He i recombination, it 
is not necessary to include a very large number of shells. Unlike 
in the case of hydrogen recombination, the exponential tail of 
He II — > He I recombination, which potentially is the most sen- 
sitive to the completeness of the atomic model, is entirely buried 
by the large number of free electrons from hydrogen. In addition, 
practically no ionized helium atoms remain after recombination, 
although in the case of hydrogen a small residual fraction re- 
mains. This is because there are significantly more electrons per 
helium atom than for hydrogen, such that freeze-out for helium 
occurs at an exponentially lower level. 

Our results suggest that for He ii ^ He i recombination the 
inclusion of 5 shells is already enough to capture the evolution of 
Xe during this epoch with precision better than 0.1 %. This pre- 
cision is sufficient if one is interested in cosmological parameter 
estimation from the angular power spectra (Cf's) of the CMB. 
However, still rather significant modifications of the recombina- 
tion history can be expected in particular from the feedback of 
He ii-photons and pro bably other physical pr ocesses that were 
omitted here (e.g. see iSwitzer & Hiratall2007al) . 



5. Bound-Bound helium recombination spectra 

In Fig. [U we present the main result of this paper, namely the 
complete bound-bound helium recombination spectrum, arising 
both during the epoch of He in -^ He ii (5000 < z < 7000), and 
Hen — » He I recombination (1600 < z < 3000). For compari- 
son, we also included the results obtained in our previo us com- 
putati ons for the H i bound-bound recombination (.Chluba et al.l 
l2007h . and added the additional line appearing as a consequence 
of the re-processing of He i photons in the continuum of hydro- 
gen, as described below (see Sect. 15. 2t . Also the 2s two-photon 
decay continua for all cases are shown. There are two important 
issues to be mentioned here: 

(i) First, the helium spectral features (both for He ii and He i) are 
significantly narrower than those of the hydrogen recombina- 
tion spectrum. This is due to the fact that for He i, recombina- 
tion occurs significantly faster due to the inclusion of the hy- 
drogen continuum opacity, and in the case of He ii, because 
its recombination occurs much more close to Saha condi- 
tions in the first place. Even the inclusion of Doppler broad- 
ening due to electron scattering, as described in Sect. 15.51 is 
unable to change this aspect. As a consequence, both recom- 
bination spectra contain clear features in the low frequency 
domain (v ~ 1 GHz), where the H i spectrum is practically 
featureless. This increase in the amplitude of variability of 
the recombinational radiation at low frequencies might help 
to detect these features in the future. 

(ii) Secondly, the He i recombination spectrum displays two neg- 
ative features, at positions v ^ 145 and 270 GHz. This is 
qualitatively different from the case of the hydrogen and He ii 
spectra, where the net bound-bound spectra appear in emis- 
sion. As we will discuss below, the reason for these features 
is directly connected with the dynamics of recombination. 
They are associated with transitions in which the lower state 
is effectively "blocked" for all downward transitions, such 
that faster channels to the 1 ' So level are provided through 
energetically higher levels. 



We now discuss in detail some particular aspects of the recom- 
bination spectra. 

5. 1 . Importance of the hydrogen continuum opacity for the 
bound-bound He i spectrum 

In Fig. |5] we show the comparison between the He i spectrum 
in three cases, namely the case with hydrogen continuum opac- 
ity assuming full redistribution of photons in the resonance; the 
case with hydrogen continuum opacity assuming partial redistri- 
bution; and the case without the inclusion of the hydrogen con- 
tinuum opacity in the computation. 

The width of the lines, which is directly connected with the 
duration of the recombination process, is significantly smaller 
when including the effect of hydrogen continuum opacity. In ad- 
dition, the peaks of the lines are shifted towards lower frequen- 
cies (i.e. higher redshifts). As a consequence, the spectrum has 
a richer structure as compared to the case of hydrogen, since the 
overlap of lines is smaller 

The full redistribution and partial redistribution spectra are 
very similar in the low frequency (v ^ 30 GHz) region. However, 
at higher frequencies, several differences appear In particular, 
for the full redistribution computation there are three negative 
features instead of two. The spectrum for the case without con- 
tinuum opacity is much smoother than the previous two, and 
presents only one negative feature. In addition, in this frequency 
regime (and specially for the strong feature at v ^ 2000 GHz) 
is seen that, due to the different speeds of the recombination 
process, the lines appear displaced towards lower frequencies 
(higher redshifts) as we move from the lower to the upper panel. 

For the high frequency region (v > 100 GHz), we present 
a more detailed direct comparison in Fig. |6] between the cases 
of partial redistribution and the one without continuum opacity. 
In this figure, a linear scale in the vertical axis is used in order 
to emphasise the existence of the negative features. One can see 
that the relative contribution of the different lines is strongly al- 
tered. In general, all emission appearing above v ~ 500 GHz is 
suppressed, while at lower frequencies, some lines are enhanced. 
We can understand these changes as follows: the contributions 
appearing at v > 500 GHz correspond to the «'Pi - I'So-series 
of neutral helium, the spin-forbidden transitions directly con- 
necting to the ground state (n^Pi - I'So), and the two-photon 
decay of the 2 'So singlet-state. The first two series contribute 
most to the strong feature at v > 2000 (see Fig.|7]for some more 
detail), while the broad two-photon continuum dominates the 
spectrum in the vicinity of v ~ 1000 GHz. When the hydrogen 
continuum opacity is included, in our current implementation of 
the problem, only the Hei 2'Pi - I'So and He: 2^1 - I'So 
intercombination-transition are directly affected, i.e. via the in- 
clusion of APesc, whereas all the other lines are modified only 
indirectly because of to the change in the recombination dynam- 
ics and the rela tive importanc e of different escape channels. 

Figure 8 of Wong & Scotti (120071) shows that without the hy- 
drogen continuum opacity the Hei 2'Pi - I'So channel defines 
the rate of recombination at z > 2400, while at z < 2400 the 
2^Pi - 1 'So spin-forbidden and, to a smaller extent, the 2'So two- 
photon decay channel dominate. They computed that of all elec- 
trons that reach the ground state of helium 39.9% go through the 
He 1 2'Pi - 1 'So transition, 42.8% pass through the 2^^Pi - 1 'So 
spin-forbidden transition and only 17.3% take the route via the 
2 'So two-photon decay channel. 

In our computations including the hydrogen continuum 
opacity we have to keep in mind that there is a fraction of elec- 
trons that reach the helium ground due to continuum absorption 
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Fig. 5. A comparison of the He i recombination spectrum for rimax - 30 with (top and middle panels) and without (bottom panel) 
the inclusion of the effect of hydrogen continuum opacity. Upper panel corresponds to the case of full redistribution of photons in 
the He I 2 ' Pi - 1 ' So -resonance, while the middle panel corresponds to our final (fudged) computation for the partial redistribution 
case (see text for details). Ordered in this way, from top to bottom we have progressively a slower recombination, so all the features 
become broader. In all three cases, a solid line indicate positive values, while a dot-dashed line indicates negative ones. Two-photon 
decay continuum 2 ' So - 1 ' So is also included as dotted line in both panels. 



by hydrogen, which then lead to the emission of additional pho- 
tons in the H i recombination spectrum. Direct integration of the 
total number of photons in the neutral helium spectrum around 
V ~ 2000 GHz yields A^^ = An/c J dvAIy/(hv) ~ 0.46A^He, while 
for the 2' So two-photon decay spectrum we have ~ 0.16A^He 
photons. Moreover, one can compute the number of photons in 
the newly generated hydrogen Lyff line (see Figure|9]below), ob- 
taining A^y(Lya') ~ 0.44A^He- These numbers show that ~ 90% of 
all electrons that reach the ground state of helium pass through 
the He I 2'Pi - I'So and Hei 2-'Pi - I'So intercombination- 
transition. The 2 'So two-photon channel only allows ~ 8% of 
all helium atoms to recombined, and ~ 2% go through the other 
n'Pi - I'So and2^Pi - I'So spin-forbidden transitions. 

We note that the modification of the dynamics of He i recom- 
bination is influencing the relative amplitude of other lines, such 
as the 3'D - 2'P (6680 A) transition, which is strongly ampli- 
fied, or the 3^D - 2^P (5877 A) transition, which now appears in 
absorption. We will discuss these transitions in detail below. 

To end this subsection, we remind again that those features 
in the vicinity of v ~ 2000 GHz arising from the Hei 2'Pi - 
I'So and Hei 2^Pi - I'So intercombination-transitions will not 
be observed in the real spectrum, because of feedback processes 
connected with HI continuum absorption at lower redshifts will 
take away these photons and will produce additional distortions 
in the HI spectrum. 



5.2. Importance of the hydrogen continuum opacity for the 
bound-bound H i spectrum 

In Fig. [8] we show how the hydrogen recombination spectrum 
is modified because of the additional free electrons appearing 
due to the absorption of He i-photons in the hydrogen contin- 
uum absorption. Most of them recombine after a very short 
time through the main channel of hydrogen recombination at 
high reds hifts, which is th e Ly-a transition (e.g. see Fig. 10 in 
Rubino-Martin et alj|2006l) . producing a "new" hydrogen Ly-a 
feature at z ~ 1870. In Fig. |9] we show the shape of that line 
separately. However, as explained in Sect. 13.31 the small addi- 
tion of electrons to the continuum also produces changes in the 
rest of the hydrogen spectrum, as shown in Figure |8] In some 
cases (see e.g. the high-frequency wing of the Paschen series) 
the changes are important at the level of 10 percent. This feature 
is due to the new Ho- feature. On average, the new bound-bound 
H I spectrum is slightly higher in amplitude, as a consequence 
of the additional photons appearing in this process. Because the 
re-processing of Hei photons occurs at high redshifts (above 
z - 1800), the two-photon continuum line is practically un- 
changed. 



5.3. Negative features in the He i spectrum. 

One of the most interesting results of our computations is the ex- 
istence of two negative features in the He i recombination spec- 
trum. In order to identify the transitions which contribute most to 
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Fig. 6. A comparison of the He i recombination spectrum for 
"max - 30 close to the He i 2'Pi - 1 'So line, with (solid curve) 
and without (dashed curve) the effect of hydrogen continuum 
opacity. The case with continuum opacity corresponds to the par- 
tial redistribution (fudged) computation. The two-photon decay 
continuum 2'So-l'Sois also included in the spectra. 



Fig. 8. A comparison of the Hi recombination spectrum for 
«max = 10 at high frequencies, with (solid curve) and without 
(dashed curve) the hydrogen continuum opacity was included in 
the treatment of the He i atom. The 2s two-photon decay contin- 
uum is also shown in both cases. 
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Fig. 7. Transitions in He i atom from the n -2 shell to the ground 
state, in particular the Hei 2'Pi - I'So and Hei 2''Pi - I'So 
intercombination-lines. This figure was obtained using our re- 
sults for the Hmax = 20 computation, including the effect of 
hydrogen continuum opacity in the partial redistribution case 
(fudged solution). 
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Fig. 9. Hydrogen Ly-a recombinational line. When including 
the effect of continuum opacity and solving simultaneously the 
evolution of hydrogen and helium recombination, most of the 
electrons which are taken from the helium by neutral hydrogen 
atoms re-appear as a hydrogen Ly-a line at v x; 1320 GHz [z ~ 
1870). 



those features, in Table[T]we provide a list with the position and 
amplitudes of all the individual lines which are found to be neg- 
ative at an amplitude smaller than -1 x 10"^^ J nT^ s"' Hz"' sr"' 
from our «max - 20 computation. To help in the discussion, we 
also provide in Table|2]a list with the position and amplitudes of 
all the positive individual lines which are found to have an am- 
plitude larger than 1 x 10"'"^ Jm"^ s"' Hz"' sr"' from the same 
computation. 

We now discuss in detail each one of these two negative 
features. Figures [TO] and [TT| present them separately, together 
with the main contributors according to the list of transitions in 
Table [1] and Table |2] For completeness, we also discuss in this 
subsection the feature which is associated to the 2'So — > 2'Pi 
singlet-singlet transition. Figure [12] presents the contribution of 



this line to the total spectrum in that region. Although this partic- 
ular line is negative, the overall spectrum in the region is positive 
due to the added contribution of other lines. 



5.3.1. First negative feature (v * 145 GHz). 

As Fig.fTOlindicates, the largest negative contribution is coming 
from one of the 10830 A fine-structure lines. The appearance 
of this feature in absorption can be understand as follows. The 
channel connecting the 2^Si triplet level with the singlet ground 
state via the two-photon decay is extremely slow (~ 4x10"^ s"' ), 
and therefore renders this transition a "bottleneck" for those 
electrons recombining through the 2^Si triplet state. Since the 
electrons in the 2^Pi triplet level can reach the I'So level via 
the much faster (~ 177 s"') Hei 2-'Pi - I'So intercombination- 
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Table 1. Positions and amplitudes of all the negative lines in the 
He I recombination spectrum for n^ax = 20 with peak intensities 
smaller than -1 x 10"^^ J m"^ s ' Hz ' sr '. We show the corre- 
sponding terms for the lower and upper states, the peak intensity 
at the minimum, the central frequency (vq) as observed today, 
the redshift (Zmin) at the minimum, and the wavelength for that 
transition in the rest frame (/Irest)- 



Fine structure Hel 10830A line (2'P-2'S) 



Lower Upper 



Aly at minimum vq 

[Jm--s-'Hz-'sr-'] [GHz] 



^min 



2'So 


2'P, 


3^D2 


4'F, 


3^D, 


4'F3 


3^02 


5'F, 


43D2 


5'F. 


23s, 


23Pi 


23p, 


3^0, 


23P, 
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-1.1 X 10"^'* 
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[A] 



78 
86 
86 
125 
40 
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273 
272 
86 
86 
40 



1855 
1875 
1875 
1875 
1875 
1905 
1870 
1875 
1870 
1875 
1875 



20590 
18690 
18690 
12790 
40380 
10830 
5877 
5877 
18690 
18690 
40380 



Table 2. Positions and amplitudes of all the positive lines in the 
He I recombination spectrum for n^ax = 20 with peak intensities 
greater than 2 x lO"^"^ Jm-^ s-' Hz-' sr-^ We show the corre- 
sponding terms for the lower and upper states, the peak intensity 
at the maximum, the central frequency vo as observed today, and 
the wavelength for that transition in the rest frame (/liest)- 
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2^8, 


23p, 
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Fig. 10. First negative feature in the He i spectrum. The largest 
negative contribution is coming from the 10830 A fine-structure 
line. Note that the upper level has fine structure, so we present 
all three possible values of J. This figure was obtained using our 
results for the n^ax - 20 computation, including the effect of 
hydrogen continuum opacity. 



He I 3P — 2D transitions 



4x10"^' 



E 1 XI 




150 200 250 300 350 400 

1/ [GHz] 

Fig. 11. Second negative feature in the Hel spectrum. The 
largest negative contribution is produced by the 3D-2P (triplet- 
triplet) transition. This figure was obtained using our results for 
the «max - 20 computation, including the effect of hydrogen 
continuum opacity. 



transition, this provides a more viable route. On the other hand 
the 2^Po and 2^P2 do not have a direct path to the singlet ground 
state. But as one can see from Table |2] this restriction can be 
avoided by taking the route 2^Po/2 -^ 2^Si -^ 2^Pi -^ I'Sq. 
The relative amplitude of these lines seen in Fig. [TOl also sug- 
gests this interpretation. 



5.3.2. Second negative feature (v * 270 GHz). 

The second overall negative feature in the bound-bound Hei 
recombination spectrum is mainly due to the superposition of 
the negative 5877 A and positive 6680 A-lines (see Fig. fTTl i. 
Here it is interesting that in Table |2] no strong positive triplet- 
singlet transition appear, which actually starts with a 3''D2-state. 
However, as Table[T]shows a strong flow from the 3^D2-state to 



higher F-level exists, again permitting electrons to pass to the 
singlet-ground level because of singlet-triplet mixing. This then 
also contributes to the close-by emission feature via the chain 
23Pi -» 3^D2 -^ 4^3 ^ 3'D2 -^ 2iPi. 



5.3.3. The spectrum in tine vicinity of v ^ 80 GHz. 

Fig. fT2l shows that in this spectral region, there is a clear low- 
intensity feature in the overall spectrum, which is produced 
by the 2' So —> 2'Pi singlet-singlet transition, that contributes 
as a negative line. Comparing the 2' So two-photon decay rate 
(^2'So-i'So)' with the transition rate to the 2'Pi, shows that at 
z ~ 2500 the latter is a factor of 2 x lO'* larger Therefore, when- 
ever escape in the He 1 2'?, - 1 ' So Hne substantially controls the 
rate of helium recombination, the 2 'So —> 2'Pi singlet-singlet 
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Fig. 12. The He i recombination spectrum in the vicinity of the 
20590 A line. This figure was obtained using our results for the Fig. 13. Influence of electron scattering on an initially narrow 
n^ax = 20 computation, including the efi^ect of hydrogen contin- Une for diff'erent emission redshifts. The vertical lines indicate 
uum opacity. the epochs of recombination at which most of photons are emit- 

ted. 



transition appears in absorption. As explained above, account- 
ing for the hydrogen continuum absorption this is the case at all 
redshifts of importance. 

However, the situation is a bit more involved, since sev- 
eral other transitions contribute to the negative (3''D2 -^ 4'F3, 
3^D3 -^ 4'F3, 3^D2 -^ 4^F3 and B^^Ds -^ 4^3) and positive 
(4IF3 -^ 3ID2 and A^Y-i -^ S'Dz) centered at v ^ 85 GHz. 
The superposition of these lines then yields an oscillatory fea- 
ture between 80 and 90 GHz, which although always positive, 
still shows the clear signature from the 20590 A line. Here 
it is important to realize that several triplet-singlet transitions 
are involved, allowing triplet atoms to decay further to the sin- 
glet ground state. This emphasises the importance singlet-triplet 
mixing for the spectrum, and in particular well-mixed levels like 
the low »F-states and beyond (mixing angle ~ 45°, see Table 
1 1.12 in [Prakeld 19961) ) provide very attractive routes. 



5.4. The He u-recombination spectrum 

The recombination history of He in is the o ne which is most 
close t o the Saha-solution (e.g. see Fig. 15 in Switzer & Hirata 
l2007bl) . Therefore the release of photons occurs during a shorter 
period than in the case of H 11 and He 11 recombination. In com- 
parison to hydrogen the release of He ill recombination pho- 
ton happens at roughly 4 times higher redshift and temperature 
(roughly 1400 for hydrogen as compared with 6000 for He iii). 

As Fig. [1] shows, the high frequency feature always appear 
on the red wing of the corresponding hydrogen lines. However 
at low frequencies the oscillatory feature drop out of phase with 
the hydrogen lines. It is also interesting to see that the He 11 and 
He I bound-bound spectra show constructive (at v > 10 GHz) 
and also destructive (v ~ 2 - 5 GHz) interference. As men- 
tioned above, this fact strongly increases the probability to ob- 
serve these features in the future. 



5.5. Changes of the spectra due to electron scattering 

The procedure to approximately include the effects of pho- 
ton scattering off free electron is outlined in the Appendix [C] 
However, here we neglect the corrections to the recombination 



history and recombination spectra arising from the changes in 
the escape of photons from the optically thick resonances, but 
these are expected to be rather small. 

In Fig. [13] we show the comparison of the Doppler broad- 
ening and recoil term for different redshift. During the epoch 
of hydrogen recombination Doppler broadening is less than 1 %, 
while it exceeds ~ 2% during He 11 — > He i recombination, and 
reaches ~ 7% at the beginning of He ill —> He 11 recombina- 
tion. We also show the strength of the recoil term for Xy - \ and 
Xy - 40. The latter case provides an estimate for the equivalent 
of the Lyman-Q' line of the corresponding atomic species. One 
can see that during hydrogen recombination the recoil term is 
completely negligible. During He 11 — » He i recombination the 
Hei 2'Pi - I'So line is shifted by < 1% and only for the Hen 
Ly-ff line the recoil shift is comparable with the broadening due 
to the Doppler term. We therefore shall neglect the recoil term 
for the hydrogen and He 11 ^ He i recombination spectrum. 

Figure [14] shows the importance of the Doppler and recoil 
term for the He in -^ He 11 bound-bound recombination spec- 
trum. At low frequencies Doppler broadening strongly lowers 
the contrast of the quasi-periodic intensity pattern, while as ex- 
pected the recoil term is not important at there. Similarly, the 
high frequency features are slightly smoothed out due to the 
Doppler effect, but the recoil term only becomes important for 
the He 11 Lyman- and Balmer-series photons. However, since in 
this work we have not yet included the re-processing of He 11 
photons in the continuum of He i and also the feedback absorp- 
tion by hydrogen, we shall not consider the corrections to the 
He ni-recombination spectrum due to the recoil term any further 
A more complete treatment of this problem will be left for some 
future work. Also, given the overall uncertainty in our model of 
the neutral helium atom we did not include the effect of Doppler 
broadening for the He n-recombination spectrum. 

6. Discussion 

In this Section we now critically discuss the results presented 
in this paper for the helium recombination spectrum. We ex- 
pect that an overall ~ 10% - 30% uncertainty is associated with 
our modelling of neutral helium, while neglected physical pro- 
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Fig. 14. Influence of electron scattering on the He ill —> He ii 
bound-bound recombination spectrum. The upper panel shows 
the changes at low frequencies, where recoil is negligible. The 
lower panel illustrates the importance of the recoil term for the 
He II Lyman- and Balmer-series. 



cesses are expected to lead to modification of the resulting he- 
lium bound-bound spectra by ~ 10% - 20%. 

We would like to mention that in addition to the aspects dis- 
cussed below another ~ 30 - 40% rather smooth contribution to 
the total recombination emission can be expected from the free- 
bound components of hydrogen (Chluba & Sunvaevll2006ah and 
helium, possibly with stronger signatures at high frequencies. 

6.1. Uncertainties in our modelling of the helium atom 

6.1 .1 . Completeness of the atomic model 

The probably largest uncertainty is connected with our model 
of neutral helium. First of all, for our final bound-bound He i 
spectrum (see Fig. [T]) we only included levels with n < 30. 
As is known f rom computation of the hydrogen recombina- 
tion spectrum (iRubino-Martin et al.l l 2006t IChluba & SunyaevI 
l2006al : IChluba et al.ll2007l) at low frequencies the level of emis- 
sion strongly depends on the completeness of the atomic model. 
Therefore we expect rather significant modifications of the He i 
spectrum at frequencies below a few GHz. Computations includ- 
ing up to 100-shells or more are probably necessary. In the case 



of the bound-bound He ii spectrum, for which we already in- 
cluded /-resolved 100-shells, the results are probably converged 
with s imilar accuracy as the one for hydrogen (see Chluba et al] 
(|2007|) for discussion). 

We also have not considered quadrupole transitions in our 
computations. Ho wever, from the typic al values of the oscillator 
strengths (see e.g. Cann & Thakka r 2002), one would expect that 
their inclusion should produce small changes. 

6.1.2. Photoionization cross-sections 

The next large uncertainty is due to the use of re-scaled hy- 
drogenic approximations for the high-n photoionization cross- 
sections. We expect diff'erences at a level of 10% - 20% due 
to these. Here in particular the exact frequency dependence of 
the cross-section may influence the importance of stimulated re- 
combinations, which become important for excited levels. Even 
for n - 5 there are notable differences, when using hydr ogenic 
formulae instead of the fits by Smits ( 1996) and Benja min et al.l 
(Il999h . Moreover, we find typical differences of the order of 5- 
10% (and in some lines 20%) between these fits and the pho- 
toionization cross-sections obtained from the TOPbase database 

(ICuntoetal.ll 19931) . 

Furthermore, as shown in Fig. 1 of IChluba et al.l (l2007l) for 
hydrogen, due to the strong dependence of the Gaunt-factor on I, 
for large n most of the recombinations actually go via levels with 
small I. In particular the S and P states of neutral helium should 
still have significant non-hydrogenic contribution for n > 10, 
which we did not account for in our model, again yielding a 
rather large uncertainty for He ii — » He i recombination. Due to 
the full hydrogenic character for the wave functions of the He ii 
atom, there is no significant uncertainty due to the cross-sections 
for He III -^ He ii recombination. 



6.1.3. Energies and transition rates 

In terms of level energies and transition rates, our model of the 
neutral helium atom is probably accurate on a level of 1 % - 10%. 
The main uncertainty is connected with t he neglect of singlet- 
triplet mixing for n > 10. As Table 11.12 in lDrakd d 19961) shows, 
for n = 10 the P and D states are still nearly orthogonal, while 
the F states are reasonably mixed, and mixing is practically com- 
plete for all the other levels. However, there are reasons why this 
may not be of so extreme importance: the highly excited levels 
(n > 10) are already very close to the continuum. Therefore the 
route via the continuum leads to a quasi-mixing of the high lev- 
els. In addition, the cascade of electrons to lower levels, where 
mixing is fully included, is very fast, such that no significant 
blocking of electrons in the higher levels is expected. However, 
the emission of low frequency photons probably will be under- 
estimated. Here, a more rigorous analysis is required. 



6.2. Additional physics missing in our computation 

6.2.1 . H I continuum opacity 

As discussed in Sect. |3] the hypothesis of complete redistribu- 
tion is not valid for the Hei 2'Pi - I'So line. This as sumption 
was usually very good in the context o f hydrogen lines JGrachevI 
119891: iRvbicki & dell' Antoniol |1994|) . in particular due to the 
presence of a huge amount of CMB blackbody photons, which 
allow electrons to pass to higher levels while they are undergo- 
ing a resonant scattering event. However, it is not the case here. 
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because of the additional continuum opacity due to small traces 
of neutral hydrogen during helium recombination. 

In order to efficiently compute the escape probability in the 
Hei 2'Pi - I'So line, we make an ansatz about its redshift de- 
pendence. This means that we fudge the no redistribution so- 
lution to the escape probability with a certain factor which is 
obtained from a diffusion code which treats in detail the escape 
problem. Although this simplification may introduce errors in 
the frequency spectrum of the order of few percent, we consider 
it acceptable given the uncertainty that we have in the atomic 
model and the photoionization cross-sections. 

We also made the simplification of assuming the va- 
lidity of complete redistribution for the Hei 2''Pi - I'So 
intercombination-line. The exact treatment of the escape of pho- 
tons in this line may lead also to differences in the spectrum 
of the o rder of ten percent. In agreement with Switzer & Hiratal 
(l2007bl) our more detailed computations (Chluba et al. 2008, in 
preparation) show that here electrons scattering plays an inter- 
esting role. 

Finally, we stress that in this paper, we only considered 
the detailed computation of the deviations of the escape prob- 
abiUty with respect to the Sobolev approximation for the He i 
2^Pi - I'So transition and Hei 2-'Pi - I'So intercombination- 
line. Inclusion of all the other n'Pi - I'So and spin-forbidden 
transitions may lead to corrections of the order of ~ 10 % aswell. 



6.2.2. Hei continuum opacity 

The absorption of He ii Lyman-a photons by the small fraction 
of neutral helium atoms during He ill —> He ii recombination 
will lead to the appearance of additional He i-photons, just like 
in the case of hydrogen (see Sect. 15.21 ). But since the number 
of photons emitted in the He ii Lyman-o- line is comparable to 
the total number of helium nuclei, this will be a notable change. 
Most obviously the He ii Lyman-a line will nearly disappear 
In addition this will accelerate He in — > He ii-recombination, 
bringing it even closer to the Saha-solution. 

6.2.3. Feedback processes 

As mentioned in Sect. |4] for He ii — > He i-recombination one 
of the probably most important processes that we neglected in 
our computations so far is feedback. As we have seen in Sect. |5] 
(e.g. Fig. |7|l, the total number of photons emitted in the He i 
2^Pi - I'So transition is comparable with those coming from 
the spin-forbidden 2^Pi — > I'So line. The former has an en- 
ergy that is larger by Av/v ~ 1%. Therefore one expects the 
Hei 2'Pi - I'So photons to interact with the Hei 2^Pi - I'So 
intercombination-resonance after a very short period of red- 
shifting. The maximum of the Hei 2'Pi - I'So line appears at 
z ~ 2550 (see Table|2]i, such that the bulk of these photons reach 
the spin-forbidden transition at Zf ~ 2520. At this redshift the 
optical depth in the spin-forbidden line is < 1, so that this feed- 
back will not be complete. Still one should check this process 
more carefully. 

As mentioned above, there is some pure continuum absorp- 
tion, far away from the resonances where resonance scattering 
can be neglected, which is not included into our program. This 
process should also lead to the re-processing of the remaining 
Hei 2'Pi - I'So and Hen Lyman-c photon, such that practi- 
cally only the hydrogen Lyman-a line will survive at the end, but 
potentially with interesting traces of the recombination history 
from earlier epochs. Also the feedback due to photons emitted 



in th e Hei n'Pi - I'Sp-series (se e lSwitzer & Hiratal ( l2007ah and 
also IChluba & SunvaevI (l2007al) for more detail) and similarly 
for He II, should lead to some modifications. However, these are 
expected to be rather small. 

6.2.4. Two-plioton decays 

The simplest addition to the two-photon processes is the 
inclusion of stimulated emission a s suggested earlier for 
hydrogen (Chluba & S unvaevI |2006b) and also included by 
Hirata & Switzer ( 2007l) for helium. These should modify the 
2s two-photon continua at the percent level. However, we have 
shown that when accounting for the effect of hydrogen contin- 
uum absorption on He ii ^ He i recombination, only 8% of all 
helium atoms reach the ground state via this channel. Hence one 
does not expect large changes in the He i recombination spec- 
trum. 

For the epoch of He in — » He ii recombination this may be a 
bit different, since electrons in higher levels will feel the change 
in the support of the levels from below, because at that time 
the two-photon decay channel is more important. In our com- 
putations ~ 44% of all electrons reach the ground-state of He ii 
via the two-photon channel, and the rest passes trough the He ii 
Lyman-series. In addition one should include the re-absorption 
of escaped heliu m Lyman-a photons by the two -photon process 
as discussed by JKholupenko & IvanchikI (|2006|) for the case of 
hydrogen. 

Also one co uld think about the two-photon decays 
from h igher levels (JSwitzer & Hiratall2007al : IChluba & SunvaevI 
l2007cl) . but both in terms of additional photons and increase of 
the overall rate of recombination one expects corrections at a 
level less than 1%. 



6.2.5. Collisional processes 

In the computations for this paper, collis ional processes have not 
been taken into account. As discussed in lChluba et alj (|2007|) for 
the case of the hydrogen recombination spectrum, because of the 
large entropy of the Universe, collisional processes only mod- 
ify the populations of the hydrogen levels for very high shells. 
In that paper it is shown that /-changing collisions need to be 
included only for shells above n ^ 30 - 40, while n-changing 
collisions can be neglected even for shells as high as « « 100. 

In the case of helium recombination, the same qualitative be- 
haviour is expected. Although in this case there are more elec- 
trons and protons per helium atom than in the case of hydrogen, 
we still expect a small effect of collisions, and which mainly 
would affect the high-« shells, i.e. it would only have an impact 
on the low frequency tail of the recombination spectrum pre- 
sented in Figure[T] A detailed consideration of the importance of 
collisions on the results will be left for a future work. 



7. Conclusion 

We have presented detailed computations of the contributions 
to the cosmological recombination spectrum due to bound- 
bound transitions in primordial helium. The re-processing of 
Hei2'Pi - I'So andHei2^Pi - I'So intercombination-photons 
by neutral hydrogen has been taken into account, yielding a sig- 
nificant acceleration of He ii — > He i and hence much more nar- 
row features than without the inclusion of this process. In ad- 
dition, some hydrogen photons are released prior to the actual 
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epoch of hydrogen recombination around z ~ 1 100 - 1500, with 
distinct traces due to the hydrogen Ly-a transition (see Fig.[TJ. 

Probably the most interesting result is the presence of two 
negative features in the He ii ^ He i recombinational spectrum. 
This is qualitatively different from any of the other spectra dis- 
cussed so far (H i and He ii). One of those negative features is 
associated to fine-structure transitions in neutral helium. 

As illustrated in Fig.lT] the total cosmological recombination 
spectrum contains non-trivial signatures of all recombination 
epochs. We emphasize this fact in Figure [15] were we present 
a detailed view, using linear intensity scale, of three regions in 
the recombination spectrum covering the low, intermediate and 
high frequency domain. Although the relative number of helium 
to hydrogen nuclei is rather small (~ 8%), constructive and de- 
structive interference of the oscillatory emission patterns at low 
frequencies, and strong non-overlapping lines at high frequen- 
cies may provide a unique opportunity to determine some of the 
key cosmological parameters, and to confront our current pic- 
ture of recombination with experimental evidence. Interestingly 
the signatures due to helium may allow a direct determination of 
its relative abundanc e, much before the first app earance of stars, 
and as pointed out in lSunvaev & Chlubal ( l2007b . these measure- 
ments do not suffer from limitations set by cosmic variance. 

As we outlined in Sect. |6] several neglected processes have 
to be studied in connection with helium recombination, in order 
to obtain definite predictions, possibly with additional revisions. 
Nevertheless, all the results presented here strongly depend on 
our understanding of atomic physics and the processes in the 
early Universe. Currently in particular the data for neutral he- 
lium may still not be sufficient. Here help from atomic physicist 
is requked in order to increase the availability of more complete 
accurate and user-friendly atomic data, in particular for the pho- 
toionization cross-sections and transition rates. 

All the numerical predictions for the recombinational 
lines obtained in this paper, which were used to produce 
all the figures in this paper, can be downloaded from 
http://www.iac.es/galeria/jalberto/recomb. 
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Appendix A: Voigt-prof ile 

Evaluations involving the well-known Voigt-piofile (e.g. see lMihalaal 19781) : 
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are usually extremely time-consuming. However, convenient approximations can 
be given in the very distant wings and also close to the center of the hne. In Eq. 
)A. U X = ^^ denotes the dimensionless frequency variable, and the Voigt- 
parameter and Doppler- width of the line are defined by 
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Fig. 15. Relative contribution of the Hi, He i and He ii bound- 
bound recombination spectra to the total spectrum at high- (top), 
intermediate- (middle) and low- frequencies (low). Helium re- 
combination spectra (both He i and He ii ) modify the shapes of 
the existing hydrogen features, shift the peaks positions and in- 
troduce new features which represent changes of 30-40% respect 
to the H I recombination spectrum alone. 



respectively. Here vg is transition frequency and A21 the Einstein coefficient for 
spontaneous emission for the considered resonance, mne ~ 4mp is the mass of 
the helium atom. Note that for the spin-forbidden 2'Pi — I'So transition the 
Voigt-parameter is ~ 170 times smaller than for the Hei 2'Pi - I'So transition. 



For \x\ < 30 we use t he approximatio n based on the Dawson integral up to 
sixth order as described in lMihalasI 119781 Sect. 9.2, p. 279). In the distant wings 
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of the line (|.v| > 30) we apply the Taylor expansion 



3 - 2a^ 15 - 2()a^ 105(1 - 2a^) 
1 + :;— + ; + :—, 



2x^ 



Ayi* 



(A.2) 



'1 — I'So and spin-forbidden 2 Pi - 1 So transition we checked 
that the Voigt function is represented with relative accuracy better than 10" in 
the whole range of frequencies and redshifts. Using Eq. lA.2l on the red side of 
the resonance one can approximate the integral^ = \ ^ i^(x')djc' by: 



Swings — 

nx 



1 + ■ 



3-2a^ 



3-4a- 



6,r 



.2 



4,r 



^ 15(1- 



Ic?) 



(A.3) 



as long as x < —30. Since a ~ 10" , this shows that the distant wings only a very 
small fraction of photons is emitted. Using the symmetry of the Voigt-profile one 
finds ^(.v) = 1 -^(-.v), such that Eq. lA.3l is also applicable for x > 30. 



Appendix B: Computation of the APesc 

In this appendix, we focus on some numerical issues which are relevant for the 
evaluation of the integral in Eq.|2] which gives the escape probability in the case 
of complete redistribution of the photons in the resonance. 



B. 1 . Analytical approximation of Afesc 

For the He i n'Pi - 1 ' So-series photons ts 3> 1 at epochs important for helium 
recombination. In particular, for these one can also expect that tl 2> t^ at the 
relevant redshifts. Therefore the integrand of Eq. J2) will cutoff exponentially 
due to the factor e"'"'- , while the term 1 - e"'"'^ does not change extremely fast. 
For the spin-forbidden transitions this condition is not fulfilled. 
Using x(.^) as variable one can rewrite the integral (2) as 



AP^, 



f dX f \se-^^'"''[l-e-^'^-'''^">]dAx', (B.l) 



with Ax' = x' ~X- The problem now is the computation of rc(^, A^'). FromEq. 
(3b) for [v' - v]/v <s; 1 it follows 



rc(x,/) » r]c(v)\y' - y] = ;7c(v)Avd[.v' - x\ 



(B.2) 



with T\c - 
wiite Ax 



\(y) 



Assuming that A.v 



- X is sufficiently small one may 



Ax' l<l>(x). It is easy to estimate that this approximation is always 
very good within the Doppler core, while it is rather crude in the distant wings. 
Inserting this into Eq. jB.l) it is possible to cany out the inner integral analyti- 
cally, yielding: 



AP" 



Jo 



d^ 1 



-Tsd-.*-) , 



Kixm- 



-[Ts+fcCr)](i-.4-)r 



(B.3) 



, where both v and x are func- 



with fcOr) = ric(v) Avo/(f>ix) and k(x) = f^TfiAF)' 
tions of ;(^. Numerically this integral is much easier to take than the full 2D- 
integral given by Eq. (2)- As we will show below this approximation works very 
well at low redshift. 



B.2. Numerical evaluation of APesc 

To carry out the 2-dimension integral (2) is a cumbersome task. We used different 
integrators from the NAqf] and CuB4j-library and only after several independent 
attempts finally reached agreement. It is extremely important to include the full 
domain of frequencies, extending the integration to the very distant wings of the 
resonance. However, due to the extreme dift'erences in the Sobolev escape prob- 
abilities of the Hei «'Pi - I'So-series (rs ~ 10^) and intercombination lines 
(ts ~ 1) different numerical schemes are required. In order to assure conver- 
gence of our numerical integrator we made sure that we can successfully repro- 
duce several limiting cases, for which analytical approximations can be found. 
In particular, we reproduced the approximation given in the previous paragraph. 
In order to understand the behaviour of the integrand in Eq. (2) we define 
the inner integral 



Fix) = TS 



£ 



>(x)cl>(x')e 



-rL(.v..v')| 



^-T.(-V.y)l 



Ax', 



(B.4) 



See http://www.nag.co.uk/numeric7l 



Download available at: |http://www.feynarts.de/cuba/| 
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Fig. B.l. Inner integrand of APesc, as defined by Eq. (IB. 4b . for 
the He I 2'Pi - I'So line (upper panel) and the spin-forbidden 
2^Pi - I'So transition (lower panel). For all redshifts, the solid 
line corresponds to the full numerical integral, while the dot- 
dashed line represents the ID approximation, as deduced from 
Eq. (Ib3]i. 



such that APesc = jL„ ^ (f) d-V. For the He i n' Pi - 1 ' So-series one always has 
Ts 3> 1, such that the exponential factor e"'"'-''' ' = e"'"sAA' is dominating the 
behaviour of the integrand. For given x we numerically determined the frequency 
y such that e"'"L(>'.>' ) < f ^ typically with e ~ 10"'*. This turned out to be rather 
time-consuming, but we found that even within the Doppler core a sufficient 
estimate for .v' could be obtained using the wing expansion of the Voigt-profile, 
yielding the condition 



X -X-- 



1 1 



Ts 



(B.5) 



However, here we typically chose e ~ 10" , in order to achieve agreement with 
the more rigorous treatment. For the spin-forbidden transitions this simplification 
is not possible, since none of the exponential factors really saturate. The full 
range of frequencies x < x' had to be considered in this case. In practise we 
never went beyond lO'' Doppler width. 

In Fig. IB. II we show F(x) for the Hei 2'Pi - I'So line and the He: 
2'Pi - 1 ' So intercombination-transition. For the He: 2'Pi - 1 ' So line the inner 
integrand becomes very broad at high redshifts, with significant contributions to 
APesc out to several thousand Doppler width, while it becomes rather narrow at 
low redshifts. However, we found that the outer integral for APesc nearly always 
has to be carried out within a very large range round the line center. One can also 
see that the approximation of F(x) following from Eq. IB. 3) works extremely 
well at low redshifts. For the Hei 2 Pi — 1 So intercombination-transition the 
main contributions to AP^sc always come from within the Doppler core and the 
wing contribution is ~ 10"^ - 10"^ times smaller For the full numerical integra- 
tion it therefore is possible to restrict the outer integral to a few hundred Doppler 
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width. In Fig. lB.ll one can again see that at low redshift the approximation from 
Eq. jB.3t works veiy well. 



Appendix C: Inclusion of line broadening due to 
electron scattering 

The photons released in the process of recombination scatter repeatedly off mov- 
ing electrons. In the low temperature limit this process can be described using 
the Kompaneets-equation. Neglecting the small difference in the photons and 
electron temperature at redshifts z SS 500 and introducing the dimensionless fre- 
quency variable Xy = hv/kTy, neglecting induced effects and the recoil term for 
an initial ly narrow line, centered at Xyp and released at Zcm, on e can find the 
solution iZeldovich & Sunvaevlll96ftlSunvaev & Titarchukll 198(1) 



AIiXy,Z ■■ 



-^bopplei- ,,3 



Xy AI(Xyfi, Zem) 



■y.o 



(ln.vy+3v(.'-lnAy q)- 

-fy.O 



(C.l) 



where A/(j:yo>Zcm) denotes the spectral distortion at frequency Xyo and redshift 
Zcm without the inclusion of electrons scattering, and the Compton y-parameter 
is given by 



y.iz) 



r 



' kT^ 



CN^CTJ 



m^c^ H{z'){l + z') 



dz' 



(C.2) 



Note that AI(Xyfj, Zi:m)/^^ q '^ An(Xyo,z), where An is the difference of the pho- 
ton occupation number from a pure blackbody, is independent of redshift. As 
Eq. iC.U show due to the Doppler etfect the line broadens by (compare also 
iPozdniakoY et all 19791) 



Ay| 

^ I Doppler 



2 V.Vc In 2. 



(C.3) 



and shifts towards higher frequencies by a factor e^^'' . 

Also including the recoil term, to our knowledge, no analytic solution to the 
Kompaneets equation has been given in the literature. However, to estimate the 
effect on the spectrum one can neglect the diffusion term and finds that the line 
shifts by 



AvI 



' -Ve JCy.O 



(C.4) 



towards lower frequencies. There is also some line broadening connected with 
the recoil effect, but it is completely negligible in comparison with the Doppler 
broadening. From Eq. iCAi is it clear that the high frequency lines will be af- 
fected most. Note that in contrast to the recoil term Doppler broadening is inde- 
pendent of the initial photon frequency. 

To account for the effect of Doppler broadening on the final spectrum one 



only has to integrate Eq. jC.lt for fixed Xy over all possible Xyo for a given 
transition. The emission redshift Zcm of the contribution can be found with vq/v = 
1 + Zcm- Afterwards in addition the sum over all transitions has to be carried out, 
yielding the final results. To estimate the influence of the recoil etfect one can 
simply add the recoil shift of each line to the frequency before summing over all 
possible transitions. 
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